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Abstract 

The relation between the covariant Euclidean free-energy and the canoni- 
cal statistical-mechanical free energy F'~' in the presence of the Killing horizons is 
studied. is determined by the covariant Euclidean effective action. The defi- 
nition of F*-" is related to the Hamiltonian which is the generator of the evolution 
along the Killing time. At arbitrary temperatures F^ acquires additional ultravio- 
let divergences because of conical singularities. The divergences of are different 
and occur since the density ^ of the energy levels of the system blows up near 
the horizon in an infrared way. We show that there are regularizations that make 
it possible to remove the infrared cutoff in After that the divergences of F^ 
become identical to the divergences of F^. The latter property turns out to be 
crucial to reconcile the covariant Euclidean and the canonical formulations of the 
theory. The method we use is new and is based on a relation between ^ and heat 
kernels on hyperbolic-like spaces. Our analysis includes spin and spin 1/2 fields 
on arbitrary backgrounds. For these fields the divergences of ^, and F^ are 
presented in the most complete form. 
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1 Introduction 



There are two approaches how to describe quantum thermal effects in the gravitational 
field. The approach by Gibbons and Hawking defines the partition function of the 

system as an Euclidean path integral. It enables one to express the free energy F^[g, (3] in 
terms of the effective action iy[5f,/5], as F^[g,(3] = l3~^W[g, P]. Functionals are 
given on Euclidean manifolds TW/j with the period (3 in the Euclidean time r. (3 is con- 
sidered as the inverse temperature. The fields are assumed to be periodic or antiperiodic 
in r, depending on their statistics. The Gibbons- Hawking approach is a straightforward 
generalization of the finite-temperature theory in the Minkowsky space-time. Its advan- 
tage is that it is manifestly covariant in the Euclidean sector and enables one to consider 
the gravitational field on the equal footing with matter fields. This approach is espe- 
cially convenient for the thermodynamics of black holes it reproduces the entropy, 
temperature and other characteristics of a black hole in the semiclassical approximation. 

When the space-time is static statistical-mechanical quantities can be also described 
in a canonical way. The canonical partition function has the form 

= Tr e"^^ , (1.1) 

where the Hamiltonian Ti is the generator of the time evolution of the system. Z*^ is well 
defined when 7i is a normally ordered operator In this case the free-energy takes the 
form 

F<=^[g, (3] = In Z = r]p-' H du^^ In (1 - r^e"^") , (1.2) 

JO duj 

where r] = +1 for Bose fields and = — 1 for Fermi fields. is the density of eigen- 



values uj of quantum-mechanical Hamiltonians of the fields 0-|T3- The advantage of 
definition (|1.2| ) is that it is given in accordance with the unitarity evolution of the system. 
However, as distinct from the Gibbons-Hawking approach, it is not manifestly covariant. 

There is no special terminology^ to distinguish and F*-^. In this paper we call 
F^ and F*^ the covariant Euclidean and the canonical free energies, respectively. The 
corresponding formulations of the finite-temperature theory will be called the covariant 
Euclidean and the canonical formulations. Sometime we will also say "Euclidean" instead 
of " covariant Euclidean" , for simplicity. 

For static space-times without horizons comparing F^ and F^ shows that these 
functionals differ only by the vacuum energy which is not included in F*". Thus, in this 
case the Euclidean and canonical formulations are, in fact, equivalent. 

In space-times with Killing horizons the quantum theory has a number of specific 
properties. On one hand, in the Euclidean formulation there is a distinguished value 
of the period f3, corresponding to the Hawking temperature, for which A^^ is a regular 

-'^ Allen called exp{—(3F^) and exp(— /SF*^) the "quantum" and "thermodynamic" partition func- 
tions. These names are not very suitable at least because the both partition functions are essentially 
quantum. 
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spaceQ. At other values of f3 the space Aip has conical singularities which result in 
additional ultraviolet divergences [15]- [0. On another hand, the canonical formulation 



runs into difficulties because the time evolution is not defined at the bifurcation surface 
of the Killing horizons. The density ^ of the energy levels in Eq. (|1.2| ) blows up near 



the horizon |jT9|-||2T| at any temperature in an infrared way. 

As a result, in the presence of horizons the Euclidean and the canonical free energies 
look different and finding the relation between them becomes a problem. This problem 
has not been analysed before and our aim is to investigate it for the case of scalar and 
spinor fields in some details. Comparison of the covariant Euclidean and the canonical 
formulations is important for different reasons. The main of them is statistical-mechanical 
interpretation of black hole thermodynamics which can be naturally defined in the frame- 
work of the Gibbons-Hawking approach 

The paper is organized as follows. In Section 2 we describe the Euclidean and the 
canonical formulations of statistical mechanics of scalar and Dirac fields on static back- 
grounds. The ultraviolet divergences appearing in because of conical singularities are 
given in Section 3 in the most complete form. We use the dimensional and Pauli-Villars 
regularization procedures. In Section 4 we develop a method how to find the divergences 
of the density of levels We first apply this method to study ^ in the presence of 
the spatial cutoff near the horizon. Then we show, in Section 5, that in dimensional and 
Pauli-Villars regularizations the spatial cutoff can be removed and one can define ^ on 
the complete space. In these regularizations the divergences of ^ have the ultraviolet 
character. It means that the corresponding divergences of the covariant free energy 
coincide exactly with the divergences of the Euclidean free energy F^. In Section 6 we 
discuss a hypothesis that (for spinors and scalars) the entire bare functionals F*" and 
F^ must coincide. We illustrate it with some examples. Technical details are given in 
Appendixes. In Appendix A we remind the reader how to relate the canonical free energy 
on ultrastatic spaces to the effective action. Appendix B is devoted to the calculation of 
the spinor heat coefficients on conical singularities, some of the coefficients represent the 
new result. 



2 Definitions and basic relations 

Let us consider scalar fields described by the Klein-Gordon equation and spinor fields 
ip described by the Dirac equation, 

{-V^V^ + ^R + m^)^ = , (7'^V^ + m)^ = , (2.1) 

■^This is the property of nonextremal black hole backgrounds. The extremal black holes will not be 
considered here. 
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where R is the scalar curvature and are the covariant derivatives, defined according 
with the spin of the fields^. The Dirac 7-matrices 7^ = (7", 7") obey the standard 
commutation relations {7'^,7'^} = 2g^^ , 7° is the anti-Hermitean matrix. It is supposed 
that the space-time is static 

ds^ = goodt^ + gahdx^dx^ , a, 6 = 1, 2, 3 . (2.2) 

The component Qqq is a nonpositive function and (700 = — 1 at spatial infinity. The 
temperature measured at infinity is On space ( [^.2D equations ( ^.11 ) can be rewritten 
in the form 

- + = , ~i^\idt + Ha)iJ = Q , (2.3) 

Hl = \go,\{-V''Va-w''Va + m^ + iR) , (2.4) 

Hd = «7o (7"(Va + ^Wa) + . (2.5) 

Here Va is the covariant derivative computed with the help of the metric Qab of the 
three-dimensional surface of constant time t = const. We denote this surface B. Index 
a is up and down with the help of gab- The vector Wa = |Valn|(7oo| is the vector of 
acceleration. The operators Hg and H^, are called the one-particle Hamiltonians because 
their eigen- values coincide with the frequencies of one-particle excitations 0. To calculate 
the canonical free energies Ff [i = s, d) with the help of Eq. ( |1.2| ) one has to know the 
densities ^ of the energy levels of Hs and Hd- 

Let us define now the Euclidean free energies for the fields described by Eq. (|2.1j ). To 
this aim we consider the Euclidean manifold Aip which is the Euclidean section of the 
Lorentzian geometry ( |2.2|) 

ds^ = g^rdr^ + gabdx^dx'' , < r < (3 , (2.6) 

where grr = \goo\- The Euclidean effective actions for the fields and ip are 

Ws[g,f3] = hogdetLs , Wd[g, (3] = -log det , (2.7) 

= -V^V^ + + , = 75(7'^V^ + m) . (2.8) 

It is assumed that Wi are regularized functionals. The operators Ls act on scalar fields on 
which are periodic in r, act on spinors which change the sign when r is increased 
by jS. The Euclidean matrix 7,- is i'-fo, the matrix 75 anticommutes with the other 7's 
and it is "normalized" as 7I = 1. Both operators ( |2.8| ) are Hermitean with respect to the 
standard inner product {(p,<f') = J ^d^x. According to Eq. ( |1.2| ), the canonical 

■^We define the spinor derivative as = 9^ + F^, where = |[7'^,7''] V^'V^V^a is the connection 
and VI are the tetrades. 

It is easy to check that iJ^ and Hd are Hermitean operators with respect to the following inner 
products (<^',(/)) = /g V(3).g|5ooh^ {<\>')*<\>, W ,i>) = /g d^x {^'Y^ where (3)^ = detg^,, and 

('0')"'" denotes Hermitean conjugation. 
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free energy vanishes at zero temperature. It is convenient to define the Euchdean free 
energy so that it have the same property, i.e., as 

F,^[g,P]=p-'WA9,P]-E^[g] , (2.9) 

E^[g]= Yim((3-^W,[g,f3]) . (2.10) 

The quantities E^\g\ have the meaning of the vacuum energy. Note that F^\g^(i\ and 
Ef\g,fi\ are covariant functional of the metric, because the Euchdean actions PV^^',/?] 
are covariant at any values of /5. 

Our aim is to find the relation between F^ and F'^ . The important property of F'^ is 
that it can be represented in the form similar to Eq. ( |2.9| ) 

Ff\g,(5\ = (5-^W\g,i5\-E%\ , (2.11) 

E%\=\mi[p-^W\g,(i^ . (2.12) 

The functionals M/j are the following effective actions 

= ^logdetL, , = -logdetLrf . (2.13) 



The operators Lj are related to Lj, Eq. ( ^.8] ), by the conformal transformations 

= e'^ , = e-i'^Lrf ei" , (2.14) 



where e~^'^ = grr- Representation (|2.11|) is well known @]-[|ri|l but for the sake of com- 
pleteness we give its derivation in Appendix A. The key point is that Li are expressed in 
terms of the one-particle Hamiltonians 

L, = El-dl , (2.15) 

Ld = 757.(^d + 9.) , Ll = Hl~dl , (2.16) 
Hl = e-''Hy , H^ = e-I^ H^e-^'^ , (2.17) 
where •jr = e.'^lr- Note that according to Eqs. (|2.17| ) the spectra of the operators Hi and 



Hi, Eqs. ( p.4|) , ( p.5|) , coincide, which means that Hi are simply an another representation 
of Hi. The operators Li act on the fields on the ultrastatic space /Ajs with the metric 

df = dr^ + gabdx^dx^ , < r < (3 (2.18) 

which is conformally related to metric (|2.6|), g^j^i, = gfiu/\goo\- One can show that 

Hd = ^loir^a + e-''m) , (2.19) 

where = '^g^" ■, and 

= -V"Va + e-^'^m^ + Vi . (2.20) 
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The derivatives Vq, are defined with the help of the metric Qab of the surface r = const 



in ultrastatic space ( |2.18|) . This surface is conformally related to B and we denote it B. 
The "potential terms" Vi are 
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(2.21) 
(2.22) 



Here = (0, w^) is the four vector of acceleration and R is the curvature of the ultrastatic 
background 



R = e"^" {R + 6{V^w^ - w^w^)) 



(2.23) 



The formulation of statistical mechanics in terms of the theory on ultrastatic spaces was 
developed by Dowker and Kennedy ||^ and by Dowker and Schofield [T^,||14|. Relations 
(g:T9D-(g:23D coincide with those used in Refs.[|13l,[ll| . 



The actions Wi and Wi are determined by the conformally related wave operators, 
see Eqs. ( p.l4| ). In static space-times without horizons the renormalized functionals Wi 



and Wi differ by the anomalous terms computed in [13|,|n]. These terms, however, are 
proportional to j3, so they result in the difference between vacuum energies ( p.lO| ) and 
( p.l2|) . The Euclidean and the canonical free energies in this case coincide. When there 
is a horizon the conformal transformation to the ultrastatic metric becomes singular and 
this case requires a special analysis. 



3 Covariant Euclidean formulation: conical singular- 
ities and divergences 

To begin with we describe the class of space-times with Killing horizons which will be 
discussed here. We assume that the metric ( |2.2| ) near the bifurcation surface S has the 
following form 

ds^ = goo{0,p)dt^ + dp^ + -fpg{e,p)dePde'' , p,g=l,2 . (3.1) 

In this representation the location of S is determined by the equation p = 0, 6^ are the 
coordinates on this surface. We also assume that near p = the components of the metric 
are decomposed as 

goo{e,p) = -k'p' (l - lp'R^J^Ae) + 0(p^)) , (3.2) 

7U(^,p) = %,i9)+p'hM + 0{p') , (3.3) 

where is the metric tensor on S. The constant k, is called the surface gravity. Metrics 
which obey the properties ( p.l|) - (|3.3|) correspond to static nonextremal black holes. It 
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can be shown that the quantities Rijij, Ru are the projections of the Riemann and Ricci 
tensors calculated on E on the directions normal to this surface. Namely, 



Ri 



R n^n" 



Rijij 



(3.4) 



Here nf are two unit orthonormal vectors orthogonal to S and the summation over the 
indexes i,j = 1, 2 is assumed. It can be also shown that 



Near S the Euclidean section M.^ oi space-time looks as 

ds^ ~ K^p^dt^ + dp^ + %qdePde'^ , < r < /3 



(3.5) 



(3.6) 



This space is regular when = jSjj^ = where the constant Pj[^ is called the Hawking 
temperature. For arbitrary (3 there are conical singularities and Aip looks as C/3 x S, 
where Cp is a cone. 

As a result of conical singularities, the Euclidean free energy is divergent even 
after subtracting from it the vacuum energy. The divergent part F^^y of can be 
calculated with the help of different regularizations. We begin with the dimensional 
regularization and consider D-dimensional space-time. It is assumed that when going to 



arbitrary dimensions the background space holds its Killing structure and equations (|3.1| )- 
( |3.3|) do not change. The difference now is that the tensors R^uXp, R^u are calculated in D 
dimensions and S is a (D — 2) dimensional surface. The divergent part of the Euclidean 
free energy is 



(47r)^/2 3/€/?2 



/ 47r2 

/im^ - [vi-^^P + V2R + V?,Ri. 



(3.7) 



where D is considered as a complex parameter and the integral is taken over the bifurcation 
surface S, /j. = ^d^^'^6. We put P = ^Rijij — Ru and introduce the constants /i and 
Pk which are listed in Table 1: 



spin 


fi 


Pi 


P2 


P3 





1 


1 

60 







1 
2 




7 

480^'^ 







Table 1. 



Here is the dimensionality of the spinor representation, = 4: for Dirac spinors and = 
2 for massless Weyl spinors. As follows from Eq. ( p. 71) , in the dimensional regularization 
F^ has a simple pole at D = 4. The dimensional regularization reproduces the divergences 
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of the logarithmical type only. For this reason it is also worth studying -^Jj^ in Pauli- 
Villars regularization which usually gives all divergent terms. The Pauli-Villars method 
is based on introduction of several, say, 5 additional fields. 2 fields with masses have 
the same statistics as the original field, while other 3 ones with masses have the wrong 
statistics, i.e., they are fermions for scalars and bosons for spinors. The latter fields give 
contribution to with the sign opposite to that of physical fields. To eliminate the 
divergences two restrictions are imposed 



m 



" + E^.'-EW = ' ^^ = 2,4 . (3. 



k 



These equations can be resolved by choosing Mi 2 = A/3/i^ + m^, M[ 2 = v^/?~-Fm^, 
M3 = + m^. The divergences in Pauli-Villars regularization can be obtained from 

Eq. (|3.7|). By adding contributions of the regulator fields with the sign corresponding to 
their statistics and taking the limit D — >■ 4, which is finite due to restriction ( |3.8| ) with 
p = 2, one finds 



/ 47r2 

bfi + a \ pi -^2^2^ + + PsRv. 



(3.9) 



a = a{m,^M) = -In m'^ -Y^ln + J2HKy , (3-10) 

k r 

b = b{m,n) = mHnm^ + -T.iK?HK'f ■ (3-11) 

k r 

The parameter /i^ plays the role of the ultraviolet cutoff. The regularization is removed 
when /i — s> 00. In this limit a ~ ln(/i^/m^) and b ^ fi^ (a, 6 > 0). Thus, in general, -P^Jj^ 
includes both logarithmic and quadratic divergences. 

The derivation of Eq. ( p.7| ) is standard. In the dimensional regularization the Schwinger- 
DeWitt proper-time representation ||2^ gives the divergent part of the effective action in 
the form 

'^L[9,P] = -lf'^e--'^j^^{Bo + sB,^ . (3.12) 

Here 5^ are the Hadamard-Minackshisundaram-DeWitt-Seeley (or heat) coefficients of 
the heat kernel asymptotic expansion]^ 

Tre-^^^j^-^^{Bo + sB, + s^B, + ...) . (3.13) 

The Laplacians A look as A = — V'^V^ + X, where is the covariant derivative defined 
according with the spin; X = {1/6 — for scalars, and X = for spinors. The relation 
between A's and the operators and L^, Eqs. (p.8| ), is 



L, = As + m' , Li = Ad + m' . (3.14) 



'We do not take into account the boundaries of space-time. 
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On singular spaces A4f^ the coefficients 5^ can be represented as the sum of two terms 



Bh 



Ak + A 



f3,k 



(3.15) 



Ak have the form of the standard coefficients defined on the regular domain of A^/j, A^^k, 
k > 1, are functionals on S which appear because of conical singularities. For scalars and 
spinors the first two coefficients have the form 



A 



^13,2 = — 



37 



37 

(7^-l)piP+(7' 



1)^ 



l){p2R + P3Ri 



where 7 = %- and the numbers /i and pk are given in Table 1. Expression 



(3.16) 

(3.17) 
I) for 



follows from Eqs. (U), (CT) and (CT),(CT)-(|3T7D. 
For scalars the coefficient Ajs^i was found by Cheeger [|^, see also Refs. 
The spinor coefficient A^^i follows from the results of Refs. 



17 , 15 



251, 261. The scalar coeffi- 



cient Aj3^2 and the general structure of the higher coefficients Af^^k were analysed in Refs. 



|T^,p3,|2^. The calculation of the spinor coefficient ^4^,2 is our new result. Its derivation 
is similar to that of Ref. ||16| but has new features related to the spin. The reader can 
find the details in Appendix B. 



4 Canonical formulation: infrared divergences 

Our aim now is to investigate the divergences of the canonical free energy and compare 
them with the results ( |3.7D ,( |0| ) found in the Euclidean formulation. We begin with 
remarks concerning specific features of quantum systems in the presence of horizons. In 
this case the one-particle oscillators of fields have a continuous spectrum of frequencies 
< u < oo. The eigenvalues of the operators H^, Eqs. (|2.19D,(PT^I1D, run down to 



uj = even for massive fields, i.e., the usual mass gap is absent. This property has the 
simple explanation. Hi are given on the space B, which is the spatial part of ultrastatic 
space ( |2.18| ) related to the original space ( |2.6| ) by the conformal transformation. On the 
ultrastatic space the location of the horizon is mapped at infinity and B turns out to be 
non-compact. At the same time, the masses rrii of the fields can be neglected near the 
horizon because they are multiplied by the vanishing factor e~'^^ = \goo\- Regarding the 
"potential" terms in Hi, see Eqs. ( ^.21| ),( p.22[ ), they are constant and negative on the 
horizon, Vg = —n"^, = — and look as tachionic masses [^,|21]]. As a result, the 
densities ^ of the eigen-values of Hi blow up near the horizon in an infrared way. 

To investigate this divergence the following method can be used. By the definition, 
the trace of the heat kernel of the operator Hf is 

-H?t _ r ^, driiiu) 2, 



Tre-^-*= / du—^e-'^' . (4.1) 
Jo du 



9 



By making use of the inverse Laplace transform one obtains 

driiiuj) 



TT 



da e*°^ Tr e 



-iaHf 



The diagonal matrix elements (a;| exp(— = exp{—Hft) 



diag 



(4.2) 



are well defined but 



the traces involve the integration over the non-compact space B and diverge. The key ob- 
servation is that the divergent parts of "^^^^ can be found by making use of the asymptotic 
properties of the traces at small t, which are very well known. 
Let us begin with massless fields in the Rindler space 



ds^ = -K^p^dt^ + dp'' + dx' + dy' 



— oo < x,y < oo , p > 



(4.3) 



The Rindler horizon is the plane IR^, p is the proper distance to the horizon. The Rindler 
space can be considered as an approximation to the black hole geometry near the horizon. 
The metric on the corresponding space B is 



df 



K '^{dp^ + dx^ + dy"^ 



(4.4) 



As can be shown, B coincides with the hyperbolic manifold IH^ having the constant 



curvature R 



29 



A review of the heat kernels of Laplace operators on such 
spaces can be found in Refs. [^,|^. It is remarkable that for the massless fields the 
heat kernels are known explicitly and their diagonal elements readQ 



diag (47rt)3/2 



diag (47rt)3/2 



1 + 



(4.5) 



The number appears in the spinor case after tracing over the spinor indexes. Obviously, 
the traces of these operators diverge at p = 0. Let us restrict the integration in the traces 
by values p > e where e is a proper distance to the horizon. Such a method is called the 
volume cutoff. With the help of Laplace transform ( |4.2| ) (see Ref. [^) one easily obtains 
from (|4.5|) the divergences of the densities of levels at e — > 



dusiuJi e) 
du) 



div 



dnd{u,e) 
doj 



A 



div 



4712^3 



(4.6) 



where A formally stands for the area of the horizon []. For scalar fields our result agrees 
with previous computations by the WKB method, see, for instance Ref. |jl9|. A general- 
ization of (|4.6| ) to massive scalars was explicitly found in Refs. ||20|1,[^|. 
Let us consider how do Eqs. 



modify when the geometry deviates from the Rindler 
form. The spaces B can be approximated by only in the limit p ^ 0. For this reason, 
the diagonal elements of the heat kernels of on B are represented by the Taylor series 

^The expression for the scalar kernel is given in Ref. the spinor kernel follows from the C, function 
which is also given there. 

^Since the Rindler horizon is a plane, only the densities of levels per unit area have a strict meaning. 
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in converging at p — to expressions ( |4.5| ). Presumably, the coefficients in these series 
should be the local functions of the curvature and so to find them it is sufficient to use 
the asymptotic form of the traces. One has 

1 



-Hft 



diag (47rt)3/2 



(4.7) 



where = 1 and aj_„ are the diagonal elements of the corresponding heat coefficients. As 
before the summation over the spinor indexes in ( [4.7] ) is assumed. With the help of Eqs. 
(g:2T])-(g:23D one finds 



asA 



-2(T 



m 



(4.8) 
(4.9) 



These expressions can be decomposed in powers of by making use of Eqs. (|3.1|) -(|3y 

1 



-Rl 



,2 



2 2 



-Ru + 0{p' 



p^R,,,, + 0{p''] 



m 



+ 0{p' 



ad,i = Td— 



iff -iff 
2 2 



2m' 



0{p' 



(4.10) 
(4.11) 
(4.12) 
(4.13) 



where R is evaluated at p = 0. For the massless fields in the Rindler space Og^i = 0, while 
0^(1,1 = fd\i^^) in agreement with Eqs. ( |4.5| ). Decompositions ( [4.12| ), ( ^.13| ) are written 
explicitly up to the terms of the order p^, other terms do not contribute to the divergence 
of the traces. 

Analogous decompositions can be found for the heat coefficients aj„ with n > 2, 
however, aj^„ vanish at small p faster than p^. To see this let us consider the operators 
Ls, L\ which are given on the ultrastatic space M.^, Eq. ( |2.18| ). Z^, are related to 
the three dimensional Hamiltonians by Eqs. ( |2.15| ),( prT^ ). Because M.p = x B the 
diagonal parts of the heat coefficients of Z^, Z^ coincide with aj^„. On the other hand, Z<j, 
Z^ are related by conformal transformation ( |2.14| ) to the covariant Euclidean operators 
Ls, L^. In this case, according to Dowker and Schofield [p!3|,p^. 



-4o- 



ai,2 + Vu Jf ) 



(4.14) 



where ai^2 are the heat coefficients corresponding to L^, L^. The currents Jf have the 
formQ 



= - 



1 

45 



^Note that our signature of the metric is different from that in Refs. p^,p^. 
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15 



(4.15) 



-^^-14 
45 



ISV^Vw + 22R'"'w, 



(4.16) 

where = w^Wy and Vw = V^Wu- The coefficients aj,2 are regular at p — ^> 0. With 
the help of Eqs. ( [4.10|) ,( ^.11|) one can check that V^Jf are regular as well. Therefore, as 
follows from ( [4.14|) , 0.12 ~ e~^°" ~ p^. The other coefficients „ with n > 2 are determined 



from ai^2 by recursion relations and so ai^n vanish near E at least as fast as aj^2- 

Our conclusion is that only ai i contribute to the divergences of the traces. From Eqs. 
( [4. 2|) , (^TTI ) , (|4. 1 2| ) , ( [4 . 1 3|) we obtain the divergent part of the densities of levels 



dusiuj, e) 



div 




6 



dndjuj, e) 
dui 









- div 




2 




1 

12 




(4.17) 



(4.18) 



Here P is the quantity defined after Eq. ( p.7| ) and / is an additional infrared cutoff 
parameter imposed at a large distance from the horizon. For massive fields / ^ mT^. 

The following remarks concerning expressions (|4.17|) and (^4.18|) are in order. As com- 
pared to the computations on the Rindler space, Eqs. ( |4.6| ), there are logarithmic correc- 
tions from the mass of the fields and non-zero curvature near the horizon. Because of the 
logarithmic terms the densities ^ have a different behavior at small frequencies for scalar 
and spinor fields. Spinor density (|4.18|) is positive at small e in the whole range of fre- 
quencies. Contrary to this, scalar density ( |4.1?| ) at u; — > is proportional to n^rn^hie^ /P 
and is negative (note that curvature corrections are negligible with respect to the mass of 
the field). Such a feature indicates that for scalars the description of the modes with low 



frequencies, the so called soft modes, may need a modification Ref. [32 



From Eqs. (pr2| ) , (|4?T7|) , ( |4. 18| ) we get the divergent parts of the canonical free energies 



^fdiv[^'/^'^] 



TT 



d,div^ 
77r2 




-i]R-m' 



In 



:P + — t:t ( -Rii - —R - m' 
96/?2 V8 12 



In 



/2 



(4.19) 



(4.20) 



5760/54' 

These equations enable one to calculate the divergences of the other characteristics of 
canonical ensembles. In particular, one can see that the statistical-mechanical entropy 
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diverges near the Killing horizon and in the leading asymptotic it is proportional to 
the area A of the horizon. For scalar fields this leading asymptotic coincides with the 
WKB results by t'Hooft [|1^] and many other authors. Equations ( [4.19|) and (|4.2(]| ) also 
follow from the high-temperature expansions obtained by Dowker and Schofield |jl3[,[|l^. 
Application of these results to our case is justified because when approaching the horizon 
the local temperature unlimitedly grows. 



5 Canonical formulation: ultraviolet divergences 

Comparison of Eq. ( p.9| ) with Eqs. ( [4.19| ), P720| ) shows that the divergences of the Eu- 
clidean and canonical free energies are expressed in terms of the similar geometrical quanti- 
ties but have different dependence on the temperature. Also the nature of the divergences 
is different: diverges in an ultraviolet way while the e-divergence in F'^ has an infrared 
origin. Finally, there is one more important difference between the regularizations of F^ 
and F'-' . The ultraviolet regularizations are usually applied to operators and functionals 
but not to the background field itself. Contrary to this, the volume cutoff regularization 
makes the space incomplete and modifies the background field essentially. 

In the presence of the horizon the densities of levels have a remarkable property. 
Namely, there are regularizations of ^ which make it possible to remove the infrared 
cutoff near the horizon and to define the densities on the complete background. As a 
result, ^ acquire new divergences which correspond exactly to the ultraviolet divergences 
of the covariant Euclidean theory. 

As the first example, let us consider the dimensional regularization. The power of the 
leading divergency in Eqs. ( [4.17D ,( ^?TB| ) is determined by the dimensionahty of the space 
B. In D-dimensional space-time the leading divergence is e^~^, if D ^ 2, and ai D > 2 
one can take the limit e — > 0. After the analytical continuation to the complex values of 
D the quantities have a pole at Z) = 4. The method how to investigate ^ near the 
pole is the following. As before we use relation ( [4.2| ). By taking into account the form of 
the operators Hf, see Eq. ( p.20[ ), we can write 



e 



diag (47rt)(^-i)/2 



{ri + bi,it + bi,2t^ + ...) , (5.1) 



bs,i = K^p'{l-^)R + 0{p^) + 0{D-4) , (5.2) 



= r^y [l + (^^Ru - \Rm -\^)+ ^(/^') + 0{D - 4) 



(5.3) 



bs,2 = 0{p'')+0{D~A,p') , 6,^2 = 2^rf'«^V + 0(p^) + 0(D-4,p^) . (5.4) 

The coefficients \n are found with the help of Eqs. (|4:^),(^l2D-( ^l4D The terms 0(D-4), 
0(Z? — 4, p2) denote additions which appear in formulas (|4.12| )- (|4.14|) when D 7^ 4. We do 
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not write these terms explicitly because they do not result in singularities when D — ^ 4. 



The integration measure in the traces is obtained from Eqs. (3.1)-(|3 



X 



D-1 



K 



E JO 



After substitution of Eqs. 
the singular part of the traces 



in Eq. ( p.lD and making use of 



Tre 



-Hit 



r(i-f 



Z)-4 



div 



Tre-^^'* 



(47r)(^-i)/2 2Kt3/2 

r(i- 



m 



P 

4^ 



(5.5) 

we can derive 
, (5.6) 



m 



D-i 



R Rj 



t 



P 

4^ 



div - (AnY^-m 2ntV^ 1 + " — ) ^ " — • (5-7) 

The F-functions appear as a result of the integration over p. It should be noted that for 
the spinors the contribution of the coefficient bd^2 cancels the pole caused by the term 
rrfK2/2 in bd,i- The divergent part of density of energy levels is obtained from (|5.6|) , (|5.7| ) 
with the help of Eq. (O) 



dns{uj, D) 
du) 

dnd{uJ,D) 
duo 



r 1- 



D 



m 



D-A 



div 



div 



(47r)^/2 

r(i-f 



(47r)^/2 



K 



R Rj 



m H 

12 



(5.^ 



(5.9) 



one 



The divergence represents a simple pole at -D = 4. Finally, from Eqs. ( |1.2|) 
can find the divergent part of the canonical free energy. We do not write this divergence 
explicitly because it is exactly the same as that of the Euclidean free energy, Eq. 
computed in the dimensional regularization 



(5.10) 



This key equality can be also established in the Pauli-Villars regularization. The regular- 
ized density of states in this method is the following quantity 



driiioj, p) _ driiioj) 



du 



du 



E 



dni{u,Mk) 
du 



E 



dniju^M^] 
du 



(5.11) 



Definition ( ^.llj ) takes into account that in the Pauli-Villars method each field is replaced 
by the "multiplet" of fields. The quantities iHii^^Mjsl^ dn,{uj^M^) ^^^^ densities of levels 
of the Pauli-Villars partners and the fields with the wrong statistics give negative con- 
tributions. The number of such fields equals the number of the fields with the correct 
statistics and the leading e-divergences in Eqs. (f4.17|) , (|4.18|) are cancelled. Regarding 
logarithmical divergences Ine^, they disappear because of constraint (p.8|) with p = 2. As 
a result, regularized densities ( |5.11|) are left finite when e — > and can be defined on the 
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complete space. When the Pauli-Villars cutoff is removed {fi ^ oo 
an ultraviolet way. The divergences can be inferred from Eqs. (|5.8|) ,(l5 
account constraints 



26 + a 



dni{u),fj.) 









div 



dnd{uj,jj) 



duj 



div 



26 + a 




^ diverge in 
I by taking into 



(5.12) 



(5.13) 



Constants a and b are defined by Eqs. ( |3.10D ,( p.ll| ) and diverge in the limit of infinite /x. 
With the help of Eqs. ( |5.12| ),( |^TTB| ) one can show that 



(5.14) 



where Fg^^[g , (3 , fi] is given by Eq. (|3.9| ). It should be noted in conclusion that Pauli- 
Villars regularization of the canonical free energy was first suggested by Demers, Lafrance 
and Myers who considered a scalar field on the Reissner-Nordstrom black hole back- 
ground. The authors used the WKB method. Although our method is different the results 
for Fg^^[g,P,fi] in this particular case coincide. 



6 Discussion 



We are interested in finding the relation between the covariant Euclidean and the canonical 
formulations of statistical mechanics on curved backgrounds with horizons. Let us discuss 
first why these formulations are equivalent for spaces without horizons. As we showed in 
Section 2, the canonical formulation is equivalent to the Euclidean theory on ultrastatic 
background, Mp, Eq. (|2.18|) , conformally related to the original space-time A^/j, Eq. 
( p.6|) . The Euclidean actions Wi are determined by the operators on see Eq. 

( p.8|) . Analogously, in canonical theory the functionals Wi are determined by operators 
Li on Mp. The classical actions corresponding to these two types of operators are 

^9,'^i]= vtLiVi^/gd'^x , Ii[g,'^i]= I '^fli^iy/^d'^x , (6.1) 

where the notation (pi is used for scalars (j) or spinors ip. a. result of Eq. ( p. 14 ), 

Ii[g,iPi\=Ii[g,iPi] (6.2) 

— — 3 

for (j) = e~'^(j) and il) = e'^'^ip. The transformation from one action to another is not 
singular and the classical theories on A^/j and J^p are equivalent. In case of massless 
spinors and massless scalars with ^ = | the operators Li and Li have the same form, 
which means that the classical theories are conformally invariant. In general case this 
invariance does not exist. However, it is still possible to introduce an auxiliary conformal 
charge in the classical actions and interpret Eq. (|6.2|) in terms of a pseudo conformal 



15 



invariance |T3[,[0]. According to a common point of view the bare quantum 

actions respect the classical symmetries. Thus, for the bare regularized functionals there 
is the same equahty as for the classical action^, 

W^[9^P]haTe = ^ii9,P]haTe ■ (6-3) 

This relation is not valid for the renormalized quantities because the conformal symmetry 
is broken by the quantum anomalies |P3],[RH]. The difference between the renormalized 



actions Wi and Wi for scalar and spinor fields was found explicitly by Dowker and Schofield 



Ti\. The anomaly, which is an integral over the Euclidean space, is proportional to 
(3 and so it contributes to the vacuum energy only. As a result, the free energies and 
Fj*-^, Eqs. (|2.9|),( prTTD , coincide before and a/ter renormalization. 

In case of the horizons M.j3 and M.j3 have the different topologies, IR^ x E and x i3, 
respectively. The transformation which relates M.p and is singular on the bifurcation 
surface S. So the classical theories are not quite equivalent. On the quantum level the 
horizons result in the divergences of Ff and F^ which cannot be eliminated by subtracting 
the vacuum energy. Moreover, the divergences of Ff and F^' have the different origins. 

Our results suggest a way how can the covariant Euclidean and the canonical formu- 
lations be reconciled. We showed that there are regularizations which are applicable to 
both FJ^ and F^ . In such regularizations the canonical free energy F^' can be defined on 
the complete background M.ji and its divergences are identical to the divergences of F^ , 
see Eqs. ( [5.10| ),( |5TT^ ). When the functionals Ff and F^ are well defined, it becomes pos- 
sible to carry out the transformation from Aij^ io M.p and to interpret it as a conformal 
symmetry. In analogy with the case without horizons, we can make a hypothesis that (at 
least for scalars and spinors) the regularized bare free energies are identical 

FF[9^^]ham = F^[9,P]},^,e ■ (6-4) 



It is assumed that the both functionals in (|6.4| ) are considered in the same regularization 



To give a strict proof of this equality may be a rather difficult problem. There are examples 
which enable one its direct check. In Ref. Cognola, Vanzo and Zerbini obtained the 
free energy of massive scalar fields in the Rindler space in an explicit form, see also Ref. 
2T[| . It can be shown that these results, rewritten in the Pauli-Villars regularization. 



confirm relation ( |6.4| ). Another direct confirmation is possible in two dimensions. Two 



dimensional massless scalar fields were analysed in Ref. and these results support our 
hypothesis as well. 

Equality ( |6.4| ) enables one to apply the methods of quantum field theory to statistical 
mechanics with the horizons. As a consequence, it justifies the ultraviolet renormalization 
of statistical-mechanical quantities Our results are also important for studying the 
statistical-mechanical foundation of the thermodynamics of black holes. It was realized 



^It is true if the regularization itself does not break classical symmetries 
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in the last years that statistical-mechanical computations in this case require an off- 
shell procedure |37|, i.e., considering a black hole at temperatures different from the 



Hawking value. From this point of view, Eq. (|6.4|) demonstrates the equivalence of 
two off-shell methods, the canonical and the conical singularity methods. However, the 
comparison of the off-shell and on-shell results goes beyond the subject of this paper. As 



non- minimally coupled scalars show [^], the off-shell and on-shell computations are not 
always equivalent. 

Finally, several remarks about the range of validity of our results are in order. We 
dealt with static nonextremal black hole backgrounds. The method described in Sections 
4, 5 is applicable to the extremal black holes as well, but rotating black holes require an 
additional analysis. For the extremal black holes the density of levels ^ has the same 
property as in nonextremal case. Namely, the Pauli-Villars and dimensional regulariza- 
tions eliminate the e- divergences. For the scalar fields the corresponding calculations were 



done in Ref. [22|. Our consideration was also restricted to the scalar and spinor fields. 



The method to calculate the quantity ^ can be generalized to the fields of other spins. 
Finding the correspondence between the canonical and the Euclidean formulations for 
these fields would be an interesting extension of this work. 

Acknowledgements: 1 am grateful to V. Frolov, Yu. Gusev, and A. Zelnikov for 
helpful discussions. This work was partially supported by the Natural Sciences and En- 
gineering Research Council of Canada. 
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A Canonical free-energy and effective action 



Here we give the details how to obtain relation ( |2.11 ). We suppose that the system has 
a discrete spectrum of frequencies. Then the canonical free energy is 



F[^[g,(3]=7^,(]-'Y.d,{uj)\n{l-r^,e- 



■I3uj 



(A.l) 



where rjg = 1, rj^ = —1 and di{u) is the degeneracy of the level uj. Equation ( |1.2| ) for 
Ff[g, (3] is obtained in the limit when intervals between the frequencies u go to zero. The 
basic identities we use are 



In (l - 
ln(l + 



-I3lu 



■^ + ln/3o^+f;in (l + 



UJ 



/3' 



k=l 



-f3uj 



■^+ln2 + f:infl + 



UJ 



k=Q 



{2k + l)H^ 



(A.2) 
(A.3) 



Note that 



UJ 



— lim — 

dz 



E 

k=l-q 



k=l-q \ 



{2k + 



+ UJ' 



E 

k=l-q 



'n^{2k + qy 



(A.4) 



where g = or 1. By using the properties of the Riemann ( function we find 

'TT\2kfY' , ^ d ^ fTT^{2k+lfY' 



d °° 
lim — 

-0 dz 



d °° 
In 3 , lim — 

k=0 



ln2 . 



(A.5) 



Thus, Eqs. ( [A. 21 ) and ( |A.3| ) are represented in the form 



Q{z\uj,P)= Yl 



k=—oo 



2n 



{k + h)] +uj' 



(A.6) 



(A.7) 



where h = and Id = |. The series (|A.7| ) converge at Rez > ^, so the functions (i{z \uj, [3) 
can be defined at z — > with the help of the analytic continuation. By taking into account 
Eqs. ( [A.l|) and ( [A.6| ) we get the canonical free energy in form 



Wi[g, P] = riil- di{uj) lim P) 
2 ~^ z^o dz 



UJ 



(A.8) 
(A.9) 



where E^[g] is the vacuum energy. The quantities IVj [(?,/?] and E^[g] diverge at large uj 
although the free energy itself is finite. To see that Waifs' 5/^] is effective action ( |2.13D one 
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can regularize this functional with the help of the ( function method. It is enough to 
make the replacement 



-0 dz 



(A.IO) 



where , according to Eq. ( |A.7| ) 



^ k=—oo 



{k + h)] +uj^' 



(A.ll) 



Ci{z\P) are the generalized ( functions of the operators Li, see Eqs. ( p.l5|) , (|2.16|) . For 
operators on manifolds without boundaries (i{z \f3) can be defined as a meromorphic 



function with simple poles at 2; = 1, 2 |30 



B Spinor heat coefficients on conical singularities 



Let us consider the heat kernel of the spinor Laplace operator 
spaces Ai/B with conical singularities. For simplicity we put k = 1. 
when f3 = 2tt. The heat kernel obeys the equation 



= -V^V^ + \R on 
so M-p are regular 



{/\d)xKp{x,x' ,s) + dsKp{x,x' ,s) =^ , Kp{x,x' = 5{x,x') 



(B.l) 



where 5{x,x') is the delta function on M.^. We first describe the heat kernel on a simple 
cone which will be required for us later. It is known since Sommerfeld that heat kernels 
on Cp can be expressed in terms of the corresponding heat kernels on the plane IR^. A 
suitable generalization of the Sommerfeld representation for integer and half-odd-integer 

we can 



spins was given by Dowker |^,|^. By making use of the results of Ref. 
represent the spinor heat kernel on Cjs in the following form 
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ir^(x(r),x'(0),s) 



2i(3 J A sin ^{z + t 



-U{z)K {x{z),x'{0),s) 



(B.2) 



Here r is the polar-angle coordinate on C^. A is the contour in the complex plane which 
has two parts. In the upper half-plane it runs from (tt — e) + zoo to (— tt + e) + zoo and in 
the lower half-plane from (— vr + e) — ioo to (tt — e) — ioo. K (x, x', s) is the spinor heat 
kernel which obeys the problem (|B.1|) on IR^. The operators Kp (x, x', s) and K (x, x', s) 
correspond to the different spin structures and have the different periodicity. The kernel 
on IR^ is unchanged when going around the origin of the polar coordinate system. (There 
is no difference between the origin of the polar coordinates and other points on the plane). 
Contrary to this, the kernel on Cp changes the sign when r is increased by j3. The covariant 
derivatives on IR^ are trivial, = d^, but the the covariant derivatives on Cp are defined 
by the polar tetrades and nontrivial, = 9^ + F^. In the basis 7^ = (cri,(T2), where 
(jfc are the Pauli matrices, the spinor connection 1-form is F = —^a^dr. The role of the 
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matrix U in relation ( [B.2| ) is to make a gauge-like transformation from the derivative on 
IR^ to that on C^g, 

U{r)d,U-\T) = V^ , f/(r) = exp Qasr) . (B.3) 

In fact, when r is real U{t) is the unitary matrix which is the spinor representation of 
the rotation on the angle r. 

To find corrections to the heat coefficients from the conical singularities we follow 
the method suggested in Ref. [|I^]. According to this method, it is sufficient to work 
in a narrow domain S of the singular surface S. The rest region of does not have 
conical singularities and the heat kernel expansion on it has a standard form. S can be 
approximated as C/j x S. So here one can relate Ki3{x,x',s) with the kernel K{x,x',s) 
on the regular space J^i3=2n by the formula analogous to Eq. ( |B.2| ). The contribution 
Trgi^r^(s) from S to the trace can be written as 

TT^Kpis) = TT^K{s) + ^ [ -1^ f^d''xTT.[U{z)Kixiz),xiO),s)] , (B.4) 

where Tr^ stands for the trace over the spinor indexes. The points x{z) and x(0) are 
connected by the integral line of the Killing field dr- The two terms in r.h.s. of ( p.4| ) 
appear when contour A is deformed to a small loop around the origin and contour A' 
which consists of two vertical curves. The effect of conical singularities is related to the 
second term in ( B.4 ). The asymptotic form of K{x,x\s) is 

„-(7^(a;,a;')/4s 

K{x, x\ s) ~ -^^-^A"\x, x') J2 an{x, x')s" , (B.5) 

where a{x, x') is the geodesic distance between points x, x' and A(x, x') is the Van Vleck 
determinant. The coefficients a„ are determined in terms of the Riemann tensor and its 
derivatives. Let p be the proper distance from the points x{z) and x(0) to S. One can 



find the following relations [16 



a\x{z),x{Q)) ^ 4p2 sin^ ^ - ^p^Ri^i, sin^ z , (B.6) 

Ai/2(x(z),a;(0)) ^ 1 + ^p2i?,, sin^^ , (B.7) 

where Ra and Rijij are defined in (|3.4|). The integration measure on E can be derived 
fromEqs. (0)-(pD 



^ ./gd'^x ~ ^ ^d^-'^e J pdpdi 
For the first coefficients in Eq. ([B.5|) one finds 



(B. 



ao(x',x) ^/ + ^i?^«A.S^\x')^x" , ai(x',x) ~ --^i?/ , (B.9) 
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where / is the unit matrix in the spinor representation and S'^'' = ;|[7'^,7'^]. The matrix 
U{z) corresponds to the rotation of a vector normal to the surface E on the angle z. So 
from Eq. (p.9| ) we get 



z / 1 z \ 

Tii [U{z)ao {x{z),x{0))] ~ r^cos- f 1 + -p^ sin^ ^^ijijj 



(B.IO) 



Here is the dimensionality of the spinor representation. Now, it follows from Eqs. 
(^)-(|B7D,(^) that 

Tr, [Uiz)Kp{x{z), xiO))] 



exp 



rd- 



(4^s)D/2 



z 

cos- 



z p 1 
l + p'{ ^Rijij + ^Rii ) sin^ - + —Rijij sin^ z - —Rs 



6 



24s" 



12 



(B.ll) 

It can be shown that approximation (p.ll| ) is sufficient to find all corrections for the first 
three heat coefficients due to conical singularities. By integrating (p.ll| ) with measure 



we obtain an integral over the surface E. Its integrand is an expression linear in 
the quantities R, Ra and Rijij with coefficients proportional to cos | sin"^"^ |, q = 1,2. 
Finally, we have to integrate this expression in the complex plane, see. Eq. ( p.4|) . This 



can be done with the help of formulas [28 



dz cos ; 



1 



iP J A' sin ^z sin^ | 
dz cos I 



iP J A' sin ^z sin^ | 



180 



(7^-l)(7y + 17) 



(B.12) 



(B.13) 



2n 



where 7 
(iJ) we find 



(or, in general case, 7 
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By making use of these relations and Eq. 



(B.14) 



where Ajs^k are given by Eqs. ( p.l6|) , (|3.17|) and Table 1 for spinors. Formula (p.l4| ) 
determines the difference between asymptotic expansions on regular and singular spaces. 
It is easy to see that Ajs^k are the corrections to the spinor heat coefficients from conical 
singularities. 
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